Vector chiral states in low-dimensional quantum spin systems 
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A class of exact spin ground states with nonzero averages of vector spin chirality, (Si x Sj ■ z), is 
presented. It is obtained by applying non-uniform 0(2) rotations of spin operators in the XY plane 
on the SU(2)-invariant AfHeck-Kennedy-Lieb-Tasaki (AKLT) states and their parent Hamiltonians. 
Excitation energies of the new ground states are studied with the use of single-mode approximation 
in one dimension for S = 1. The excitation gap remains robust. Construction of chiral AKLT states 
is shown to be possible in higher dimensions. We also present a general idea to produce vector 
chirality-condensed ground states as non-uniform 0(2) rotations of the non-chiral parent states. 
Dzyaloshinskii-Moriya interaction is shown to imply non-zero spin chirality. 
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Introduction: Vector spin chirality, defined as the 
projection onto the axis of rotation (here given by z) of 
the average of the outer product of two adjacent spins, 



S J X S 7 



measures the sense of rota- 



tion of the magnetic moments in a spiral magnet. Being 
even under time reversal and odd under the inversion of 
i and j sites, this chirahty plays an important role in the 
recent study of spin-polarization coupling in multiferroic 
materials where the local dipolc moment shares the same 
symmetry properties as k^- [1]. A linear coupling be- 
tween the two order parameters is a generic phenomenon 
in spiral magnets. 

In a quite different context, an interesting observation 
was made by Hikihara et al. for the S = 1 spin chain 
with both nearest (Ji) and next-nearest (J2) neighbour 
interactions [2]. For J2I J\ larger than a critical ratio, 
the ground state was shown to possess long-range or- 
der in the chirality correlation function, (Ki^i+iKj^+i), 
as |i — j| 00. Such a novel phase, in the context 
of multiferroicity, would result in a strong coupling to 
the local dipole moment even in the absence of mag- 
netic order. A similar possibility of a non-magnetic, yet 
chirality-ordered phase was explored in the Ginzburg- 
Landau treatment of anisotropic spin models [3]. In both 
instances, the key is to reduce the symmetry of the 
Hamiltonian away from SU(2) and introduce frustration 
to suppress magnetic ordering. 

In this paper, we discuss a simple route to produce 
the vector chiral ground state in low-dimensional spin 
systems. First we look to models of Hcisenberg spin ex- 
change together with the Dzyaloshinskii-Moriya (DM) in- 
teraction ~ Si X Sj ■ z. In the one-dimensional case, one 
can consider the following Hamiltonian: 



(1) 

We choose i ^ j + l and Sf = {Sf ± iS\)l^, for a 
periodic lattice of length iV. The strength of the DM 



interaction for each (ij) bond is given by J2sin0.y. One 
can implement a site-dependent unitary rotation of the 
spins, St U^S+U} = Sfe-^^\ with the angle 6, cho- 
sen to meet the condition e*'-^*"^^-' = e*^'j [4]. A simple 
way to choose Oi is to start with 9i = 0, then choose all 
successive angles as 62 = O21, 0^ — 62 = O32, etc. accord- 
ing to 9i = J22<j<i ^j,3-i- periodic structure 
we need to require O^j^i (= 9i), which equals the sum of 
all the bond angles J2i<i<n ^i.i-ii be an integer multiple 
of 27r: 



E 

l<i<JV 



0i,i_i/27r = n{= integer). 



(2) 



Once this condition is met, it is always possible to 
"gauge away" the phase angles in the model given 
in Eq. (1) to reduce it to the XXZ Hamiltonian: 
U[{9,}]H[{9,}]U^[{0,}] = 7?xxz where U[{0,}] = ll^U^. 
The eigenstates of Eq. (1), |{6'y}), have a one-to-one 
correspondence with those of the XXZ model, denoted 
|XXZ), and given explicitly by |{%}) = W [{0i}]\XX.Z) . 

Symmetry consideration dictates that {S~ Sj~) for 
the eigenstates of Hxxz be equal to Xij, where 
Xij is a real-valued number. It then follows that 
{{^ij}\SiS+\{9,j}) = Xije''^'>^-''^'> for the eigenstates 
ll^ij}). The imaginary part of this average is nothing but 
the spin chirality, {SfSj—S^Sp, given by Xij siii{9i — 9j). 
This simple argument proves that the DM interaction in- 
duces non-zero spin chirality in the quantum eigenstates. 

For S = 1/2, through Jordan- Wigner transformation, 
the Hamiltonian (1) is mapped to a model of spinless 
fermions coupled to the gauge flux 6ij. While a per- 
sistent current will exist for general values of the flux 
4> — the criteria given in Eq. (2) corre- 

sponds to having an integer multiple of the flux quantum 
threading the ring, for which we would expect vanish- 
ing fermion current. Here, however, one must note that 
the spin chirality maps onto i{fi' fj — fj' fi) , whereas the 
gauge-invariant definition of the fermion current will be 
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i{e^^^^ fj — e^'^^^' fi). This latter quantity vanishes 
when the flux is an integer multiple of 27r but the spin 
chirality, given by i{fi' fj—fj' fi) in the fermion language, 
remains nonzero even for the integer flux case. In turn, 
{e'^'^S+Sr -e-'^'^S^Sp vanishes for |{%}) when an 
integer flux threads the ring. Our proof remains valid for 
arbitrary spin S (for which no Jordan- Wigner transfor- 
mation exists) and J2 / Ji ratio. 

The whole class of Hamiltonians given by Eq. (1) 
obeys an identical set of energy spectra regardless of 
the choice of bond angles {%}, as long as Eq. (2) 
is obeyed. We have checked this for a 4-site model 
with arbitrary {^12, ^^23, ^^34, ^41}, under the constraint 
6*12 + 6*23 + 6*34 + 6*41 = 2tt X integer, for both 5=1/2 
and S = 1 cases. The spin-spin correlation functions also 
behave in the manner predicted by the gauge argument. 

Reversing the argument, one can generate states 
of non-zero and non-uniform chirality beginning 
with the XXZ Hamiltonian by introducing a site- 
dependent phase angle 0i and rotating each spin 
accordingly: 5^ ^ UjS^Ui = Sfe^^K The 
XXZ Hamiltonian undergoing the unitary rotation 
C/t[{^J]ifxxzt/[{^i}] = H[{ei}\ becomes Eq. (1) with 
6ij = Oi — 9j. To obtain the uniform DM interaction 
one can use 9i = 9 x i where i is the local coordinate 
and require that 9N (A^=number of lattice sites) be an 
integer multiple of 27r. For the staggered DM interaction 
one can choose 9i = and 9 for even and odd sites, 
respectively. The net flux, given by Eq. (2), will be 
always zero for even A'', regardless of 9. The eigenstates, 
obtained as unitary rotations of those of the XXZ 
Hamiltonian, will have non-zero spin chirality. 

One-dimensional chiral AKLT state: We have 
presented an argument how a quantum state with non- 
zero spin chirality can be generated. The same idea 
can be applied to the well-known Affleck-Kcnncdy- 
Lieb-Tasaki (AKLT) ground states of spins for one 
dimension [5] . The discussion is most conveniently carried 
out in the Schwinger boson language where the spin op- 
erators are represented by = a|&i/-\/2, S~ = blai/\/2 
andSf = (ajoj— 6j6j)/2. Spin rotation in the XY plane is 
implemented through aj a|e*^'/^, and bl 6|e~*^'/^. 
Under this rotation, the AKLT ground state, which is 
built up of a product of bond singlet operators ^4^ = 
aj^t — fttflt, is replaced by 



i{^i})=n4[%]io)' (3) 

where = e'^^^/^alb] - e^^^-Z^fefa], and % = 

9i — 0j . The quantization rule shown in Eq. (2) is 
satisfied. The AKLT ground state can be written in 
the matrix product form[6], |AKLT) = Tr(n,5i) with 
a 2 X 2 matrix gi, and we can write down a simi- 



lar matrix product ground state for non-zero chirality, 
\{9i})=Tv{Ui9im), using 



9^\9^, 



t^2 



-iff. 



-a]b] 



(4) 



As will be shown shortly, the states given in Eq. (3) show 
non-zero spin chirality, and may be christened the "chiral 
AKLT states". The AKLT Hamiltonian undergoes the 
unitary rotation accordingly. Expressed in the Schwinger 
boson language, = Hij[9ij], the pair-wise 

Hamiltonian Hij[9ij] is given as 



H^,m - ^(6-4[%]A.,[%])(4-4,..,j..,,,.yj 
Alj[9ij]AM=2[l-S:S^ 



(5) 



The phase rotation produces e^^'^ S^S~ + e~''^^^ S~ in 
Eq. (5). Taking all 9ij = gives back the usual AKLT 
Hamiltonian. 

Having obtained the chiral extension of the one- 
dimensional AKLT state, we consider some of its ground 

state properties and the excitation energies using the 
single-mode approximation (SMA)[7]. For the ground 
state in Eq. (3), the average of S~Sj' is obtained 

from {{9,}\SrS+\{9,}) = e*^- (AKLT 1 5^5+] AKLT) = 
(l/3)e''''' {Si ■ Sj)o. Here the subscript refers to the av- 
erage with respect to the AKLT ground state. The chiral 
moment in \{9i}) follows as Kij = (1/3)(S, ■ Sj)o sin^^ = 
— (4/9) sin for nearest neighbours. The spin-spin cor- 
relation function is straightforward to work out [7] : 



m}\sf^^sf^^m) = (1/3)(S, • S,-)ocos%, 

m}\s^s^m}) = {i/3){Si.sj)o, 

{Si ■ Sj)o = 2Sij + 4(1 - ^i,)(-l/3)l'-J'l. (6) 

The identity, {{9,}\Sf S^m}) = {{9,}\Sf S^m}), is en- 
sured by the global U(l) symmetry of the chiral Hamil- 
tonian, Eq. (5). The exponential decay in the spin-spin 
correlation persists for chiral AKLT states. The ensu- 
ing SMA calculation, as well as the general argument for 
the invariance of the energy spectra given earlier, con- 
firms that the excitation gap persists for non-zero chiral 
angles. The ground state is thus non-magnetic, gapped, 
and possesses non-zero chiral moments. 

The structure factor for the uniform chiral AKLT state 
1^), where 9ij = 9 x {i — j), can be easily worked out. 
Denoting the structure factor in the AKLT state as 
s{k) = {S^S^)o = 2(1 - cosfc)/(5 -I- 3cosA;), k = -k, 
we have 



I ^ k 



\9) = [s{k+9) + s{k-9)]/2, (7) 
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E^(k,e) 




while the z-component of the string order is given 
by the usual one: O^j = S^exp[iTTj2j<k<i^k]Sj- 
The averages 



•ij< 

of the string operators for \{9i}) 

(4/9)cos% and ({^?J|Of,|{^?J) 



-4/9. 

factor cos 9ij 



3 ■ 
is 



In particular for the uniform chiral phase, the 
cos[9{i — j)] in the string order reflects 
the extra pitch angle 9 due to the helical spin structure 
introduced by the DM interaction. 

Higher-dimensional generalization: Construction 

of chiral AKLT ground states and the associated parent 
Hamiltonians are possible in higher dimensions: 



FIG. 1: SMA energies for S'^ excitations for the uniformly 
chiral AKLT state 16*) in ID chain and 5* = 1 as given in Eq. 
(8). Plots are shown for < fe < 27r and < 6* < tt. 



and {9\S^Sl.\9) = s{k). The average energy of the excited 
state 5^161) (a = x,y,z) is given by {9\SfH[{9i}]Sf\9), 
which is equal to cos 9ij {SfHSf)o for a = x, y, and H is 
the AKLT Hamiltonian. For S^, AKLT expressions are 
obtained. The excitation energies in the SMA are given 
by 



E^{k,9) = Ey{k,9) = 



f{k+0) + f{k-9) 
s{k+9) + s{k-9) ' 



(8) 



using /(fc) = (10/27)(1 — cosfc). The excitation spec- 
tra are displayed for 9 ranging from to tt in Fig. 1. 
The SMA energies for S^\9) and 5*^1 6') possess symme- 
try midcr k ^ — fc, while those for S'^\9) and S'jTl^) will 
be given by f{k ± 9)/s{k ± 6*), respectively, exphcitly 
breaking the chiral symmetry. The SMA results are also 
in accord with the general argument that all the energy 
eigenstates remain in one-to-one correspondence through 
the rotation. 

A string order parameter [8] characterizes the inher- 
ent antiferromagnetic spin-spin correlation in the AKLT 
state better than the spin-spin correlation function it- 
self, which has an exponential fall-off with the separation. 
The string order parameter in the chiral state is given by 



ix)=n(4M''i°^ 



{io) 



2S 



H'' = Y. E ^^>o- (10) 

{ij) J=2S-M+1 



Here M = 2S/ z is determined by the value of the spin S 
and the lattice coordination number z. Each (ij) bond 
carries a bond angle 9ij . With more bond variables than 
can be generated by the set of site angles, the gauge 
rotation argument of the one dimension does not readily 
apply in higher dimensions. An alternative proof is given 
as follows. 



The projector to the 



angular momentiun- J subspace[5, 

in terms 



71 in the AKLT Hamiltonian is constructed 
of the bond spin operator J^^- = (Sj-FS^O^ = 2S'(5'+ 1) -t- 
2S, • Sj. The replacement S^Sr + SrSf e'^'iSfSj + 
e-***-^ S'~S'+ in j| produces P-^ [%] in Eq. (10). To prove 
that Ix) is indeed the zero-energy ground state of in 
Eq. (10), we will show that each projector Vfj[9ij] acting 
on Ix) produces zero. 

First write 



'PiA(^ij] = VjViP^jV,^vJ, 



(11) 



ij 



where Vi is the U{1) rotation ViS^V^ = S, 



We 



S^'eKY>[i'K (cos can choose and freely as long as their difference is 



j<k<i 



equal to 9ij. Focusing on a given bond {ij), the chiral 
(9) ground state can be written out 



Ix) = E E ■ • • i4r'{blr (e^'^^4&] - e-''^^a]bl)'' {a]r^{b]r • • • |0) (12) 



4 



where the terms on the far left and far right stem from the product of j4j^^[0j,g] with only one end of {pq) connected 
to either i or j. The sum vm + tii and ruj + Uj are constrained to equal 25 — M in above. Applying the projector 
V^jlOij] on Ix) and using relation (11) we obtain 



n^-[%]ix>= E 



E 



M 



(air 



|0). 



(13) 



J — nojjYj stems from the gauge transformation V^Vj ■ ■ - ViVj applied on 



in Eq. (13) have an expansion in terms of 

— I uii tj^LL-^ is less than or equal to Jmax = 25* — M . Hence 
M , Eq. (13) is zero. Since the whole argument works for each bond (ij) we have 



Here the phase 4>j) = {rii — mi)(j)i + {rij — mj)4> 
the terms shown in Eq. (12). The state shown inside the bracket 
states for which the total momentum J,,, = S, + S, on the bond 
for projectors VijlOij] with J > 2S 
proven that -ff-^lx) = 0. 

Although the proof holds for any bond angle configuration {Oij} and for any dimension of the lattice, the gauge 
transformation introduced in Eq. (11) is only a local one, without the possibility to define the global imitary 
operator constructed as the product of Vi's. Hence, it is generally not correct to associate sin 9ij with the local 
average of the chirality except when we can decompose the bond angle as the difference of the local angles, 9ij = Oi — 9j . 



Discussion: In conclusion, we have identified a sim- 
ple and straightforward way to produce ground states of 
spins carrying non-zero vector spin chirality. The key 
idea is to start with a spin Hamiltonian whose ground 
state is non-chiral, and introduce non-uniform phase 
twists of and S~ , but not of 5?. The difference of the 
twist angle for nearby sites 6ij = 6i — 6j defines the degree 
of local vector chirality. The Dzyaloshinskii-Moriya inter- 
action also emerges in a natural way, after implementing 
the non- uniform 0(2) rotations on the Hamiltonian with- 
out the DM interaction. A simple argument shows that 
the ground state in the presence of the DM interaction 
will generally possess non-zero vector chiral moments. 

The well-known AKLT ground states of spins can be 
generalized in this way, in both one and higher dimen- 
sions. The ground state correlation properties for one- 
dimensional chiral AKLT states, in particular, can be 
readily calculated as chiral rotations of the known cor- 
relations for the non-chiral AKLT state. The excitation 
energies for the uniformly chiral AKLT state is calculated 
within the SMA and possess the gap which does not close 
as the chiral angle is varied. We have in addition iden- 
tified the string order parameter appropriate for the lin- 
ear chiral AKLT chain. Since the states we constructed 
in this paper possess nonzero chiral moment, their long- 
range ordering follows automatically. Construction of a 
different kind of chiral state, without the chiral moment 
but only long-range order in its correlations [2], will be an 
interesting challenge for the future. 

For the experiments, insulating systems having the DM 
interaction in addition to the Heisenberg superexchange, 
such as the parent compoimd La2Cu04[9]; arc the likely 
places to find ground states with non-zero vector spin 
chirality. While the lattice deformation responsible for 
the presence of DM interaction can be measured in the 



X-ray scattering, a direct, simultaneous measurement of 
the spin chirality in the same compound using the 
polarized neutron scattering [3] will highlight the corre- 
lation between the two phenomena involving the lattice 
and the spin. 
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